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Inferred cosmic-ray spectrum from Fermi- LAT 7-ray observations of the Earth’s limb 
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Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA 
reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual 
softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to 
that of protons. These newly-discovered features may offer a clue to the origin of high-energy CRs. 

We use the Fermi Large Area Telescope observations of the 7-ray emission from the Earth’s limb for 
an indirect measurement of the local spectrum of CR protons in the energy range ~ 90 GeV-6 TeV 
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(derived from a photon energy range 15 GeV-1 TeV). Our analysis shows that single power law 
and broken power law spectra fit the data equally well and yield a proton spectrum with index 
2.68 ± 0.04 and 2.61 ± 0.08 above ~ 200 GeV, respectively. 

PACS numbers: 96.50.sb, 95.85.Ry, 98.70.Sa 


Introduction. The spectrum of CRs has offered few 
clues to its origin so far. The generally accepted fea- 
tures are at very- high and ultra-high energies (see, e.g., 
Figure 1 in [1]): the so-called “knee” at a few thousand 
TeV [2, 3], the second “knee” at ~ 10 6 TeV, the “ankle” 
at higher energies [4], and a spectral steepening above 
10 8 TeV [5, 6]. It is believed that CRs below the second 
knee are Galactic, while extragalactic CRs dominate at 
higher energies [7, 8]. 

The data recently collected by three experiments, 
ATIC-2 [9, 10], CREAM [11, 12], and PAMELA [13], in- 
dicate a new feature at relatively low energy: a break (or 
hardening) of CR proton and He spectra at ~ 240 GV in 
rigidity. PAMELA claims to detect the break at 95% con- 
fidence level for both species. Below the break, PAMELA 
data agree very well with the earlier data from AMS-01 
[14] and BESS [15]. Above the break, ATIC-2 results 
agree well with those of CREAM, smoothly connecting 
to highest energy points from PAMELA. The change in 
the spectral indices for both protons and He is ~ 0.2. 

However, the break itself is observed only by PAMELA 
near its high-energy limit. Much evidence of this newly 
discovered break or flattening comes from a combination 
of data by several different experiments, which may be 
subject to cross-calibration errors. A verification of this 
new feature requires an independent confirmation, prefer- 
ably with a single instrument. Meanwhile, recent prelim- 
inary AMS-02 results * 1 do not show any feature in the 
proton and He spectra up to ~ 2 TeV and also seem to 
contradict ATIC-2 and CREAM results. In this paper, 
we demonstrate that such a measurement can also be 
done indirectly through observation of the CR-induced 
7-ray emission from the Earth’s atmosphere. 

Atmospheric 7-ray emission is mainly the result of 
hadronic CR cascades: CRs entering the atmosphere near 
grazing incidence produce showers that develop in the 
forward direction, resulting in a very bright 7-ray sig- 
nal from the Earth’s limb as seen from orbit. The 7-ray 
spectrum from CR interactions at the very top of the 
atmosphere depends only on the inclusive 7-ray produc- 
tion cross section and the spectrum of CR particles. If 
the cross section is known, the shape of the local CR spec- 
trum can be recovered from the 7-ray spectrum. How- 
ever, this method can only measure the total spectrum 
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of CRs. To deduce the spectrum of protons, the most 
abundant component of CRs, one has to assume a spec- 
trum of He, the second most abundant component. The 
contribution of the latter to the total 7-ray emission is 
~ 10-20%, depending on the energy. Therefore, accurate 
modeling of the contribution from He interactions is not 
very critical, and heavier nuclei can be neglected. 

Observations of Galactic diffuse 7-ray emission (GDE) 
have provided valuable information about CR spectra 
in distant locations [16-20] and in the local interstellar 
medium [21, 22]. Similarly, observations of the Earth’s 
limb 7-ray emission can be used to deduce the CR spec- 
trum near the Earth. In contrast with the GDE, the 
contribution from the inverse Compton scattering of CR 
electrons to the Earth’s limb emission is negligible. Fur- 
thermore, viewed from low-Earth orbit, the limb is orders 
of magnitude brighter than the GDE. Using the Earth’s 
emission is, therefore, a simpler way to derive the spec- 
trum of CR nucleons than using the GDE. 

The 7-ray emission from the Earth’s limb was first ob- 
served by the SAS-2 [23] and EGRET [24] instruments, 
but these observations were limited in statistics and an- 
gular resolution. Fermi - LAT made the first measurement 
of the Earth’s limb 7-ray emission above 10 GeV [25] and 
was able to resolve the limb profile to discriminate be- 
tween the thin and thick target regimes, demonstrating 
its capability for indirect measurements of the CR spec- 
trum. In this paper we report on the analysis of 5 years 
of Earth’s limb observations with the Fermi- LAT. 

Data and Analysis Method. Fermi was launched in 
June 2008 and spent the first few weeks calibrating the 
instruments during the Launch and Early Operations 
(L&EO) period, during which the LAT was tracking a 
few well-known bright sources, allowing the Earth’s limb 
to frequently enter the field of view (FoV) . In September 
2008, three hours of limb-stare observations were per- 
formed. This is the data set used in [25] and part of the 
data set in the analysis presented here. The additional 
part of the data set is described as follows. 

The spacecraft operates mainly in survey mode, keep- 
ing the Earth’s limb far from its boresight as it is a back- 
ground for other analyses. However, for a small fraction 
of the operating time, the LAT performs pointed obser- 
vations by following a celestial target while it is not oc- 
culted by the Earth, including while it is near the limb. 
We select this pointed data set by accepting events when 
the magnitude of the rocking angle 2 is > 52°, 2° greater 
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than that for the normal survey mode, up to August 
8, 2013. This rocking angle selection rejects the survey 
mode data, for which the Earth’s limb photons have large 
(> 62°) incidence angle. 

We avoid the geomagnetic and solar modulation of lo- 
cal CRs near the Earth by considering only 7 rays above 
15 GeV because they must be produced by CR protons 
with energies of at least (but mostly much greater than) 
15 GeV. The resulting number of Earth’s limb photons 
above 15 GeV (N>i 5 Gev) and the average incidence an- 
gles measured from the LAT’s boresight ((0lat)) for 
these datasets are in Table I. 

The data are analyzed here in the local nadir coordi- 
nates, in which ©Nadir is the angle measured from the 
nadir direction at the location of the LAT. At the LAT’s 
altitude of ~ 565 km, the physical limb of the Earth is at 
©Nadir ~ 66.7°. However, the peak of the 7-ray emission 
above 15 GeV is at ©Nadir ~ 68.1°, due to the height 
of the atmosphere and the effects of 7-ray absorption as 
discussed in detail in [25]. At ©Nadir = 68.4°, the in- 
tegrated column density for grazing-incidence particles 
is ~ 3 g cm -2 (see Figure 5 in [25]). From this angle 
outwards, the atmosphere is in the thin-target regime 
with photons produced from a single interaction, the ab- 
sorption effects are negligible, and the resulting 7-ray 
spectrum is determined by the local spectrum of CRs. 
Thus, ©Nadir = 68.4° is the inner edge for the studies 
presented here. The outer edge of the Earth’s limb is 
chosen as ©Nadir = 70.0° because the emission from ce- 
lestial sources starts to dominate for larger ©Nadir angles. 

We use the P7REP reprocessed data (see [26] and 
FSSC page 3 for details) with the P7REP_S0URCE event 
selection and the associated P7REP_S0URCE_V15 instru- 
ment response functions. We apply two additional cuts: 
©lat < 70° (reduced incidence angle) to avoid the edge 
of the FoV, which is prone to systematic uncertainties, 
and 68.4° < ©Nadir < 70° (thin-target regime) to select 
photons from the Earth’s limb. 

The background is estimated from a ring surrounding 
the Earth’s limb (80° < ©Nadir < 90°). The ring immedi- 
ately surrounding the limb was not used in order to avoid 
spill-over photons from the limb due to the LAT’s point- 
spread function (PSF). The background level shown in 
Figure 1 is small, ranging from ~ 3% at 15 GeV to ~ 5% 
at 500 GeV of the bright limb emission. 

For the P7REP data used here, dedicated simulations 
and flight data comparisons have been performed to vali- 
date the LAT responses up to 1 TeV by the LAT Collabo- 
ration. Based on these studies, we adopt an effective area 
(A e ff) uncertainty of 5% at 10 GeV, increasing linearly 


the rocking angle. 
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with the logarithm of energy to 15% at 1 TeV 4 . 

We simulate many realizations of the energy depen- 
dence of the A e g which obey the above estimated un- 
certainty to observe the effect of instrumental systematic 
error. Specifically, to obtain one realization, we gener- 
ate 3 random numbers at 10, 100, and 1000 GeV from 
Gaussian distributions for which the mean is 0 and a 
is the value of the uncertainty estimation at the three 
energy points. Cubic spline interpolation between these 
points describes the deviation of A c g from the central 
value, which would then distort the measured Earth’s 
limb spectrum in a way consistent with the systematic 
uncertainties, allowing us to evaluate the propagated un- 
certainties of the final results. This algorithm to simulate 
A e g uncertainties assumes uncorrelated errors for two en- 
ergy bins that are sufficiently far apart (larger than half a 
decade in energy) , but the interpolation results in highly 
correlated errors between nearby energy bins (see Section 
5.6.2 in [27]). 

We also correct for the angular resolution effects on 
the limb spectrum itself. The dominant effect is contam- 
ination by limb photons from ©Nadir < 68.4°, where the 
emission is brighter. The other is the leakage of pho- 
tons from the limb to each side of the defined boundary. 
Above 15 GeV, where the LAT’s PSF is narrow (68% 
containment at ~ 0.2°) and not strongly energy depen- 
dent (see Figure 57 in [27]), these corrections combined 
decrease the measured intensities by ~ 35%, depending 
on energy. The effect on the spectral index is relatively 
small compared to that from A e g. 

We determined that +2%/-5% uncertainty of the ab- 
solute energy scale (described in Section 7.3.4 in [27]) 
translates into < 10% effects on the absolute normaliza- 
tion of the spectrum, which does not alter the results 
presented here. 

To infer the CR proton spectrum from the 7-ray mea- 
surement, we use two pp — > 7 interaction models, one 
by Kamae et al. [28] ( Kamae model) and the other by 
Kachelriefi and Ostapchenko [29] ( K&O model). For 
each model, we calculate the 7-ray spectrum by integrat- 
ing a model of the proton spectrum from ~ 0.5 GeV to 
~ 500 TeV in kinetic energy. 

In our study, we assume that the atmosphere consists 
of 100% Nitrogen. This does not affect our results be- 
cause studies of proton-nucleus interactions at high ener- 
gies (e.g., [30, 31]) show that the pA cross section can be 
scaled from the pp cross section by applying an energy- 
independent scaling factor ocA 0 - 7 , where A is the atomic 
number of the nucleus. The precise scaling factor for the 
atmospheric composition is not important for this analy- 
sis because it changes only the normalization of the fitted 
proton spectrum. 


4 http: / /fermi. gsfc.nasa.gov/ssc/data/analysis/LAT _caveats.html 
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Observation type 

Start date 

End date 

Livetime (days) 

{©LAT 

’) N>15 GeV 

L&EO 

Jul 15, 2008 

Jul 30, 2008 

9 

44 ° 

967 

Limb-stare 

Sep 29, 2008 

Sep 29, 2008 

0.125 

31° 

18 

Pointed (multiple) Aug 21, 2008 Aug 8 , 2013 

90 

45° 

6762 


TABLE I. Observation types and durations for this analysis. 


Observations of the 7-ray emission from the Earth’s 
limb cannot discriminate between contributions of CR 
protons and heavier nuclei. Thus, we must rely on the 
direct measurements of the CR composition. The He 
fraction in CRs is about 6-10% by number, depending 
on energy, in our energy range of interest, so its contri- 
bution has to be taken into account, while the contribu- 
tion of the heavier nuclei can be safely neglected. There 
are a number of empirical parameterizations for nucleus- 
nucleus meson multiplicity (e.g., Appendix A in [30] and 
Eq. (3b) in [31]). These formulas give similar values for 
the ratio of aN to pN cross sections, <JaN/o'pN ~ 1.6. 
We use this number to scale the pp- interaction models 
and to calculate the relative contribution of He nuclei 
to the limb 7-ray emission. Since the contribution of 7 
rays produced by He is ~ 10-20% depending on energy, 
and most of the emission is produced by protons, the 
aN scaling uncertainty has little influence on the final fit 
results. 

To determine the He spectrum, we fit the combined 
PAMELA ([13]), CREAM ([12]), and ATIC-2 ([10]) He 
data above 50 GeV/n (102 GV) with spectral forms de- 
scribed below Equation 1. According to PAMELA mea- 
surements [13], the He/p ratio at ~ 90 GeV is 6.2%. We 
use this value together with the cross section scaling to fix 
the contribution of He to photon production at 15 GeV 
to 9.5%. We then forward-fold by varying the parameters 
of the input proton spectrum so that the resulting 7-ray 
spectrum calculated from the pp-interaction models pro- 
vides the best fit to the Earth’s limb measurement. In the 
fitting procedure, the normal Poisson likelihood function 
is maximized: 

N 

L = l[P Poisson (n° hs ,nr d ), ( 1 ) 

i= 1 

where N is the number of energy bins. Pp 0 isson is the 
Poisson probability of observing n° bs counts given that 
the model predicts n/ lod counts (7-ray model flux x ex- 
posure) for the i th energy bin. We use two models for the 
local CR proton and He spectra in the fitting procedure: 

• SPL: Single power law in rigidity. This model as- 
sumes a single power law for CR protons and He. 
For He, the index 2.73 ± 0.01 is our best-fit value 
of the combined PAMELA, CREAM, and ATIC-2 
data. For protons, we fit both the normalization 


and index to our measurement of the 7-ray spec- 
trum from the Earth’s limb. 

• BPL: Broken power law in rigidity. This model 
assumes a broken power law for both proton and He 
spectra. As before, the He spectrum is fixed to the 
best fit of the combined direct measurements, for 
which the spectral index changes from 2.82±0.07 to 
2.55 ±0.02 at 247 ±44 GV. We then fit the indices, 
break energy (Ebreak), and normalization for the 
proton spectrum. 

We evaluate the statistical uncertainties of the fit re- 
sults by fitting a large number of simulated realizations of 
photon counts generated with a Poisson distribution for 
which the expected value is the measured count in each 
energy bin. Likewise, simultaneous simulations of pho- 
ton counts and ranges of A e g, as previously described, 
give the total (combined systematic and statistical) er- 
rors. We also add in quadrature the 5% absolute energy 
scale uncertainty to the errors of the fitted energy pa- 
rameter. 

Results. The measured 7-ray thin-target limb spec- 
trum, the background-sky flux (which has already been 
subtracted from the limb spectrum), and the best- fit 7- 
ray models are shown in Figure 1. The 7-ray emission 
from the local H I 5 is scaled to approximately match that 
from the Earth’s limb and shown for comparison. As ex- 
pected, the two agree well above ~ 10 GeV, since they 
are produced from hadronic interactions by the same lo- 
cal population of CRs in the thin-target regime. Below 
~ 10 GeV, the spectra differ due to the geomagnetic and 
solar modulations of CRs in the vicinity of the Earth, 
reducing the number of CRs interacting with the Earth’s 
atmosphere. For this reason we limited our study to limb 
7 rays above 15 GeV. The approximate proton-to-7-ray 
energy conversion factor for the power-law spectrum of 
protons is 0.17 [32], The energy range of the inferred 
proton spectrum is thus ~ 90 GeV-6 TeV. 

Using the K&O and Kamae models, we obtain the re- 
sults shown in Table II. The log likelihood for the best-fit 
BPL is ~ 0.9 better than that for the best-fit SPL. To ac- 
count for Aeff systematic uncertainties, we apply Monte 
Carlo simulations to translate this likelihood ratio into 
a significance. By assuming that the best-fit SPL is the 
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FIG. 1. The 7 -ray energy spectra multiplied by E 2 ' 75 . The 
background-sky spectrum (triangles) has been subtracted 
from the limb spectrum (circles). The scaled local Hi 
emission 5 (crosses) is shown for comparison. Best-fit 7 -ray 
results from two CR models based on the K&O model [29] 
are also plotted as dotted and dashed lines. Statistical and 
total (quadrature sum of statistical and systematic) errors are 
shown as bars and bands, respectively. 

Kamae K&O 

SPL index 2.67 ± 0.05 (0.03) 2.68 ± 0.04 (0.03) 

BPL index 1 2.81 ± 0.10 (0.03) 2.81 ± 0.11 (0.04) 

BPL index 2 2.60 ± 0.08 (0.05) 2.61 ± 0.08 (0.06) 

BPL £break 276 ± 64 (55) GeV 302 ± 96 (62) GeV 
BPL vs SPL 1.0 a 1.0 a 

TABLE II. Fit results from Kamae [28] and K&O [29] models, 
shown as value ± total error ( statistical error). 

true underlying flux model, we produce ~ 2000 simu- 
lations of A e Q from the estimated errors as previously 
discussed, for each of which we generate ~ 1000 simula- 
tions of SPL realizations. We then fit the distribution of 
the log likelihood differences between SPL and BPL for 
these ~ 2M total simulations with a Gaussian function 
and evaluate how likely it is that the best-fit BPL we ob- 
tain from the actual measurement would give a log likeli- 
hood difference of 0.9 or above as compared to the SPL. 
We find that it corresponds to a significance of 1.0 er. 

We performed several cross checks to test the stability 
and consistency of the results. We studied the effects of 
using the event selection with more stringent rejection 
of residual CRs (P7REP_CLEAN), tighter incidence angle 
(©lat) cuts, and reasonable variations of the fitted en- 
ergy ranges (up to 20 GeV lower bound and down to 
120 GeV upper bound in 7-ray energy). All of these 
cases yield consistent results. 

Figure 2 shows the resulting best- fit SPL and BPL 
derived from the K&O model in comparison with di- 
rect measurements, assuming an effective atmospheric 
column density of ~ 1.0 g cm -2 , as described below. 

In order to determine the absolute normalization of 
the inferred proton spectrum, we use the NRLMSISE- 


00 atmospheric model [33] to calculate the average line- 
of-sight column density, weighted by 7-ray intensity, in 
the range studied here (©Nadir = 68.40°-70.00°) to be 
1.2 g cm -2 . Due to the exponential change of the at- 
mospheric density with ©Nadir, the evaluated density is 
extremely sensitive to the lower bound of the ©Nadir 
range. We thus empirically adjust the absolute nor- 
malization of our inferred proton spectrum to approxi- 
mately match that of direct measurements as shown in 
Figure 2 by changing the atmospheric column density 
from 1.2 g cm -2 to 1.0 g cm -2 . This is equivalent to in- 
creasing the lower bound of ©Nadir from 68.40° to 68.42° 
when we calculate the atmospheric column density. The 
small change in the effective lower bound of ©Nadir by 
~ 0.02° has many potential justifications, such as the 
LAT altitude variations which smear the precise calcula- 
tion of the target density, the atmospheric model uncer- 
tainties, and other absolute normalization uncertainties 
as previously discussed. Since our primary interest is in 
the spectral indices, the difference in normalization is of 
no importance. 

Discussion and Conclusion. Our LAT analysis, which 
employs a different technique from direct measurements, 
shows that the CR proton spectrum between ~ 90 GeV- 
6 TeV can be described equally well (~ 1 a) with the SPL 
and BPL models. The best-fit spectral indices (2.68±0.04 
for SPL and 2.61 ± 0.08 above ~ 200 GeV for BPL) are 
consistent with each other. 

We note that our best-fit SPL index is ~ 3 a from the 
value (2.801±0.007) reported by PAMELA for a lower en- 
ergy range (29-79 GeV). However, our best- fit SPL index 
for ~ 90 GeV-6 TeV is in good agreement with the fitted 
index for ~ 230 GeV-1 TeV reported by PAMELA [13] 
and with the measurements at higher energies by ATIC-2 
[10] and CREAM [12]. While Fermi-LAT results cannot 
confirm or disprove the existence of the spectral break 
itself yet, they do indicate a flatter proton spectrum at 
high energies, consistent with direct measurements by 
ATIC-2 and CREAM. 

This result is the first indirect measurement of the pro- 
ton spectrum in the energy range ~ 90 GeV-6 TeV using 
observations of the 7-ray emission of the Earth’s limb. 
Continuing observations with Ferrm-LAT will allow us 
to improve the precision of the measurement of the CR 
spectrum and extend the energy range. 
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Swedish Research Council and the National Space Board 
in Sweden. Additional support from INAF in Italy and 
CNES in France for science analysis during the opera- 
tions phase is gratefully acknowledged. I.V.M. acknowl- 
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FIG. 2. Best-fit single power law or SPL (top), and best-fit broken power law or BPL (bottom) spectrum for the local CR 
proton spectrum (solid red lines) as derived from the Earth’s limb 7-ray data using the K&O model [29] for pp-interactions. The 
total (combined statistical and systematic, neglecting errors in absolute normalization) uncertainties are the dashed red lines. 
Other direct measurements ([10], [12], [13], [14], [15]) are shown for comparison. The gray band is PAMELA’S total uncertainty. 
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